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ABSTRACT 

Based on sensitive CO measurements from HERACLES and H I data from THINGS, we show that 
the azimuthally averaged radial distribution of the neutral gas surface density (Shi+ Eh~2) in 33 
nearby spiral galaxies exhibits a well-constrained universal exponential distribution beyond 0.2 x r 2 5 
(inside of which the scatter is large) with less than a factor of two scatter out to two optical radii 
i25- Scaling the radius to T25 and the total gas surface density to the surface density at the transition 
radius, i.e., where Sjji and £h2 are equal, as well as removing galaxies that are interacting with their 
environment, yields a tightly constrained exponential fit with average scale length 0.61±0.06i"25. In 
this case, the scatter reduces to less than 40% across the optical disks (and remains below a factor 
of two at larger radii). We show that the tight exponential distribution of neutral gas implies that 
the total neutral gas mass of nearby disk galaxies depends primarily on the size of the stellar disk 
(influenced to some degree by the great variability of £h~2 inside 0.2 x i"2s). The derived prescription 
predicts the total gas mass in our sub-sample of 17 non-interacting disk galaxies to within a factor of 
two. Given the short timescale over which star formation depletes the H 2 content of these galaxies 
and the large range of T25 in our sample, there appears to be some mechanism leading to these largely 
self-similar radial gas distributions in nearby disk galaxies. 

Subject headings: Galaxies:ISM — Galaxies:Evolution 



1. INTRODUCTION 

Scaling relations are important because they are often 
thought to be manifestations of some underlying funda- 
mental physical processes. They also make it possible to 
determine properties of a system with incomplete mea- 
surements to within the scatter of the scaling. The to- 
tal neutral gas content of galaxies has long defied such 
characterizations because of the wide diversity of previ- 
ously observed gas surface density distributions (see, e.g., 
IHelfer et all 120031 This was at least partly due to the 
lack of comprehensive, resolved, extended and deep ob- 
servations of atomic (H I) and molecular gas (H 2 ) across 
the disks of a significant sample of galaxies. This limi- 
tation has been overcome by the availability of new ra- 
dio surveys providing sensitive H I and CO (the stan- 
dard tracer for H2) radio data, tracing the gas content 
of many nearby galaxies out to large galactocentric radii 
(THI NGS, HERACLES, IWalter et/ all [20081 : ILerov et al.l 
120091) . 

With these data sets in hand, we can now re-assess 
the azimuthally averaged gas distributions (radial pro- 
files) in nearby spiral galaxies. Both HI and CO pro- 
files tend to show quite distinct behavior: the HI surface 
density distributions are often depressed in the centers, 
increasing in the inner disk and slowly declining in the 
outer regions. Furthermore, the surface density in the 
outer regions can be affected by environmental factors 
such as interaction with companions or truncation due 
to the r am pressure strippin g by the hot gas in galaxy 
clusters (jCayatte et al. 1994]). As for the molecular gas, 
it tends to fall off much more sharply than the HI in 
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an approximately exponential fashion tha t varies from 
galaxy to galaxy (lYoung fc Scoviilel Il982t iRegan et al.l 
l200ltlLe7o7et^ et al.ll2008lL 

We show in this paper, however, that when the totality 
of the neutral (i.e., atomic and molecular) gas in normal 
spirals is considered based on the data from these new 
surveys, and the profiles are properly scaled, the gas ex- 
hibits a radial profile shape that is remarkably constant 
from galaxy to galaxy. The gas surface density is found 
to vary by no more than about a factor of two from this 
average profile across the optical disk and even beyond. 
We show that the scatter reduces even further when we 
eliminate galaxies that are interacting with their envi- 
ronment. 

2. DATA & METHODOLOGY 

In this paper we use radial profiles (i.e., azimuthal av- 
erages in tilted rings) of atomic hydrogen, Shi, molecular 
hydrogen (as traced by CO emission), Sh2, and the to- 
tal gas surface density E gas =EHi+5]H2- We derive the 
H 2 profiles out to one and the H I and total gas profiles 
from the centers out to two optical radii r25, which is 
defined as the 25magarcsec~ 2 B-band isophote. We will 
refer to the regime within r25 as the optical disk. 

The VLA THINGS surv ey (The HI Nearby Galaxy 
Survey, IWalter et al.l [2008) , provides sensitive (cthi ~ 
0.5 M Q pc -2 at 30" resolution), extended (~ 0.5° field- 
of-view) HI data. These data, together with new and 
archival data for some of the target galaxies, allow us 
to track the HI distribution out to 2 x r 2 5 with good 
sensitivity. 

The IRAM 30m survey HERACLES provides the dis- 
tribution of 12 CO(2-l) emission in our sa mple galax- 
ies (HERA CO Line Extragalactic Survey, ILerov et al.1 
[2001) . HERACLES provides sensitive (ct H 2 « 5 Mq pc" 2 
for the most distant, face-on systems) data covering the 
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entire optical disk out to 125- The quoted H2 surface den- 
sities in this paper include a Galactic CO-to-H2 conver- 
sion factor Xqo = 2 x 10 20 cm~ 2 (K km s -1 ) " 1 and a CO 
line ratio 7co(2 - l)/J C o(l - 0) = 0.7 fsee iLerov et al.1 
[20T1 . Because HI can be observ ed out to much larger 
radii for individual lines-of-sight, iSchruba et al.1 ()2011[ ) 
used the velocity information from the H I as a prior to 
stack individual CO spectra. This technique allows one 
to probe down to t7H2 ~ 1 M© pc~ 2 and thus to de- 
rive radial CO profiles out to the edge of the optical disk 
(~ r 25)- This represents a significant improvement over 
previous studies and we adopt their radial H2 profiles 
derived from this approach. 

We thus focus our analysis on the da ta and galaxy sam- 
ple presented in ISchruba et al.l (|2011l ). which consists of 
33 nearby, star-forming disk galaxies (excluding edge- 
on, low- mass and low-metallicity systems). All quoted 
surface densities in this paper include the contribution 
from helium (a factor of 1.36) an d are corrected for in - 
clination. We refer the reader to ISchruba et al.l f|201 lh 
for further details and references regarding sample selec- 
tion and properties, H I and CO datasets, the stacking 
technique for the CO data and the conversion of line in- 
tensities to surface densities. 

We stress that it is the combination of such deep, wide- 
field CO data, the stacking technique as well as the sen- 
sitive HI data, which permits to derive unprecedentedly 
accurate radial distributions of the atomic, molecular and 
total gas across galaxy disks. These radial profiles are 
shown for each galaxy i n the galaxy atlas in the Ap- 
pendix of Sc hruba et al.l (|2011l ) (note, however, that we 
extend their H I profiles out to 2xr2s). 

We also include radial gas profiles for the Milky 
Way for c ompari s on. W e adopt the H I and H2 pro- 
files from iDamd (|1993f Fl. adjust the latter to match 
our adopted value for Xqo and include the contribu- 
tion from helium. We assume r25 ~ 16 kpc, based 
on the observation that H II regions, planetary nebulae 
and carbon star s seem to disappear at about 2 x in 
the Milky Way dFich fc Blitd ll984t iSchneider fc Terzianl 
[19831: ISchechter et aLIHQSlT 

3. RESULTS 
3.1. Unsealed Radial Profiles 

FigureQ]shows the H I, H 2 and total gas surface density 
profiles of the 33 sample galaxies and the Milky Way. 
Radially slightly less extended versions of these profiles 
as well as the H I and CO maps they are d erived from are 
shown for each of our sample galaxies in ISchruba et al.1 
(|2011[ ). We refer the reader to the Appendix of that 
paper to match individual galaxies to particular profiles. 

The top panel of Figure Q] shows the H I profiles. 
They are generally relatively flat and show a great 
deal of scatter with little apparent regularity except 
for a concentration just below ~10M Q pc -2 for radii 
r < 10 — 12 kpc. This is the maximum HI sur- 
face density (for near solar-metallicity galaxies) be- 
fore the inter stellar medium (ISM) becomes primar- 
ily molecular dMartin fc Kennicuttl 120011; iWong fc Bht3 
120021 ILerov et al.l l2008: Bigi el et al.ll2008h . Note that the 



3 We use their H I profile based on the Kulkarni et al. ( 1982, 
KBH) rotation curve and the H2 profile following Digcl ( 1 99ll ) at 
large radii. 



profiles fill almost all of the "phase space" between 1 and 
10 M Q pc -2 at all radii. Some profiles are truncated and 
decrease steeply: these are often interacting systems or 
cluster members. In the inner parts, many galaxies show 
a more or less pronounced depression. These central "H I 
holes" are usually regions where the ISM is predomi- 
nantly molecular (the Sh2/Shi ratio incre ases roughly 
expo n entially with decr easing radius; e.g., ILerov et al.l 
120081: IBigiel et al.l f2008h . The Milky Way HI profile 
agrees well with the overall trend: H I-deficient in the 
center and relatively flat throughout at £hi ~ 5 — 7 
M Q pc" 2 . 

Almost all H2 profiles in the middle panel show an 
exponential-like decline. The scale length, however, is 
different from galaxy-to-galaxy, i.e. some profiles are 
steeper than others. Like the H I profiles, the H2 pro- 
files fill almost all of the "phase space" up to the most 
molecule-rich galaxy. In their inner parts, many galaxies 
deviate from this exponential trend and th e H2 surface 
densi ty rises much more steeply (see also iRegan et al.l 
2001). Even when averaged over many square kpc, some 
galaxies reach surface densities of many 100 M pc~ 2 in 
the inner parts, quit e similar to what is observed in star- 
burst galaxies (e.g., Kennicuttl 119981) . The Milky Way 
profile also declines roughly exponentially from r w 5 kpc 
on outward, but shows a depression at smaller radii. 
Where it declines exponentially, it does so with a scale 
length typical of many of the other galaxies. 

The plot of the total gas profiles is given in the bottom 
panel and shows that the combination of the relatively 
flat H I profiles (over large parts of the radius range con- 
sidered here) and the exponential H2 profiles leads to 
somewhat more coherent total gas profiles for many of 
the galaxies than do either the H I or H2 profiles alone. 
The total gas profiles cluster at surface densities between 
a few and a few times 10 M pc -2 between r rs 2—12 kpc. 
The Milky Way profile falls quite in the middle of the 
bulk of the profiles and thus its averaged gas distribu- 
tion matches that of most of the other galaxies well. 

We also show mean data values (red points) and asso- 
ciated la RMS scatter of the profiles in 2 kpc wide bins 
(note that the means are plotted in the middle of these 
bins) . We exclude the central 2 kpc for the binning where 
the profiles rise much more steeply than further out. The 
relative scatter ranges from ~ 60% to ~ 80% within the 
optical disks, i.e., for r < T25, and remains within a fac- 
tor of two even at large radii where the profiles begin to 
flare significantly. 

The red line shows an ordinary-least-squares (OLS) fit 
to the means and yields an average scale length of 6.1 kpc 
with standard deviation 0.3 kpc. Even though the fit is 
relatively well constrained (reflected by the small stan- 
dard deviation) , we estimate the impact of several other 
factors on the fit result by doing the following: 1) we 
add noise to each data point (mean) and re-fit, 2) we 
bootstrap (sampling profiles from the ensemble allowing 
repetition), 3) we vary the radius range to fit over and 4) 
we vary the bin size we use to derive the means. The first 
method is an uncertainty estimate based on the error on 
the mean in each bin, the second addresses how much in- 
dividual profiles drive the fit and the last two approaches 
probe the impact of our specific choices for radius range 
and bin size. Adding noise as well as varying radius 
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Fig. 1. — Radial profiles of the surface densities of atomic gas (Shi, top panel), molecular gas (Sh2i middle panel) and total gas 
(S gas =SHl+SH2i bottom panel) for 33 nearby, star- forming disk galaxies and the Milky Way. While the H I profiles are generally flat 
and concentrate just below 10 Mq pc~ 2 , the H2 profiles decline approximately exponentially. An exponential fit to the binned means (red 
points; error bars show the RMS scatter) characterizes the bulk of the total gas profiles well and the scatter about the fit remains within 
a factor of two even at large radii. 
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range and bin size lead to an uncertainty <0.5 kpc. Boot- 
strapping leads to an uncertainty of ~ 0.7 kpc, which we 
adopt as our uncertainty on the fit so that we quote the 
mean scale length of the unsealed total gas profiles as 
6.1 ± 0.7 kpc. We will show in the following that the 
agreement among these profiles can be improved further 
by scaling both axes the right way and by focusing on 
non-interacting galaxies. 

3.2. Scaled Radial Profiles 

By measuring the radius in physical units (kpc), in- 
trinsic differences in the size of the galaxies will lead to 
scatter in the profile distribution. This can be accounted 
for if the x-axis (radius) is normalized by the optical ra- 
dius (r 2 5) of the galaxy. This scaling reduces the scatter 
to 40 — 60% across the optical disk (beyond which it 
reaches up to 100% at large radii, similar to the unsealed 
profiles). This is notably less than the 60 — 80% in the 
u nsealed case above. 

Blitz fc Rosolowskvl (|2004D noted, that there is a nat- 
ural scale to the gas in normal spirals if the conversion 
from H I to H2 is governed primarily by the midplane hy- 
drostatic pressure in galactic disks. In that case, the loca- 
tion where the gas in the disk goes from being primarily 
molecular in the inner regions, to where it becomes pri- 
marily atomic in the outer regions, should occur at a con- 
stant stellar surface density. They showed that for the 28 
galaxies they analyzed, this constancy is go od to about 
40% a nd has a value of ~120 M© pc~ 2 . iLerov et all 
(2008) reached the same conclusion (although they de- 
rived a stellar surface density of 81 ± 25 M Q pc~ 2 ), us- 
ing a combination of HERACLES, THING S and BIMA 
SON G (BIMA Survey of Nearby Galaxies, iHelfer et al.l 
120031 ) data for 23 nearby galaxies. We determine this 
transition radius (i.e., where Ehi = Sjj 2 ) directly for 
each galaxy from the radial profiles (typically the tran- 
sition occurs at E gas w 14 M Q pc~ 2 ). 

In Figure [2] we show the profiles one obtains if both 
scalings are applied: the radius is normalized to the 
size of the optical disk (x-axis) and the gas surface 
density is scaled to the value of £ gas at the transi- 
tion radius (y-axis). For six galaxies the transition ra- 
dius cannot be readily determined: the ISM in these 
galaxies is either entirely dominated by H I (NGC 337, 
NGC925, NGC 4214, NGC 4559) or entirely H 2 domi- 
nated (NGC 4569, NGC 2798). We will not include these 
galaxies in this plot and the following. 

The red points in Figure [5] show mean values and la 
scatter in 0.2xr25-wide bins. The scatter increases from 
~ 35% to ~ 80% at large radii where the profiles begin 
to flare. Across the optical disks the scatter is between 
35% and 50%. While the normalization of the x-axis 
has a bigger impact reducing the scatter, the additional 
scaling of the surface density reduces the scatter further 
to 50% or less within the optical disks and to less than a 
factor of two even at large radii. The minimum scatter of 
35% is almost a factor of two improvement compared to 
the minimum scatter of 60% for the unsealed case. This 
matches the visual impression that the scaled profiles in 
Figure [5] show a remarkably tight distribution. 

The red line shows the OLS fit to the means and yields 
an average scale length of 0.48 ± 0.04r25, where the error 
is the standard deviation. Following the approach de- 



scribed in Section [3. 1[ we estimate the uncertainty of the 
fitted scale length in several ways. We find that boot- 
strapping and varying the radius range dominate the un- 
certainties and we quote our final fit result as 0.48 ± 0.06 
r25- 

The gas distribution, particularly in the outer parts, 
can be significantly affected by interactions with other 
galaxies in groups or clusters or by ram pressure strip- 
ping from a hot intracluster medium. Oftentimes these 
interactions result in truncated H I profiles, driving the 
scatter in the profile ensemble. 

In addition to the six galaxies with uncertain transition 
radii above, we further eliminate 10 galaxies that exhibit 
evidence for interactions with their environments. These 
are galaxies that show signs of strong tidal interactions 
in their HI distributions and those which are members 
of the Virgo or Coma clusters subject to ram-pressure 
stripping. These galaxies are: NGC 2146, NGC 2976, 
NGC 3077, NGC 3627, NGC 4254, NGC 4321, NGC 4579, 
NGC 4725, NGC 5194, NGC 5713 (note that NGC 2798 
and NGC 4569 would be part of this list too, but are al- 
ready excluded based on the unclear transition radius in 
these galaxies). 

Figure [3] shows the plot that results if both axes are 
scaled as in Figure [51 but after removing the interact- 
ing galaxies. For the remaining 17 galaxies, the scatter 
within the optical disks reduces to 25 — 40% and remains 
between ~ 80 — 90% beyond r 2 5. The Milky Way pro- 
file is a good match to the overall trend. The ensemble 
constitutes a tight, well-defined exponential distribution 
across the optical disk with deviations in the centers and 
flaring in the outer disks. The fit yields a scale length 
of 0.61 ± 0.03 r 2 5 and we quote our final result based on 
bootstrapping as the dominant source of uncertainty as 
0.61 ±0.06 r 25 . 

4. DISCUSSION 

Figure [3] shows that galaxies not interacting with their 
environment exhibit a tight exponential distribution of 
their total gas content when the radius is scaled to the 
optical radius and the surface density of the gas is scaled 
to the surface density of the gas at the transition radius. 
The scaling relation for disk galaxies we derive is: 

^ = 2.1xe- U5xr / r!S , (1) 

^trans 

where Etrans is the surface density of the gas at the tran- 
sition radius. The factors 2.1 and —1.65 are obtained 
from the y-intercept and slope of the exponential fit in 
Figure 3. This implies that for these disk galaxies the 
total mass of neutral gas, M gas , is given by 

M gas = 2tt x 2.1 x E trans x r 25 2 x X, (2) 

where the factor X depends on how far out from the 
center one integrates Equation [TJ Integrating out to the 
optical radius r 25 yields X = 0.18, out to 2xr 25 yields 
X = 0.31 and integrating to infinity yields X = 0.37. 

Equation [5] shows that for any given maximum radius, 
M gas depends only on r 25 2 and S tra ns- Straus does not 
vary a great deal from galaxy to galaxy and has a typical 
value of about 14 pc~ 2 (also compare ILerov et al.l 
120081: iBigiel et al.l [2008V Thus, M gas depends primarily 
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Fig. 2. — Scaled radial profiles for the galaxy ensemble: the radius (x-axis) is normalized to the size of the optical disk and S gas (y-axis) 
is scaled to the value at the transition radius, i.e., where the ISM transitions from being primarily molecular to being primarily atomic. The 
combined distribution of profiles is remarkably tight, including the Milky Way profile which is right in the middle of the galaxy ensemble. 
The exponential fit to the means is well constrained and the scatter is reduced significantly compared to the unsealed profiles in the bottom 
panel of Figure [l] 
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Fig. 3. — Scaled radial profiles as in Figure [2] but after removing galaxies showing signs of interaction with their environment. The 
distribution of these disk galaxy profiles is extremely tight, with a scatter across the optical disk of only 25 — 40% and a scale length of the 
exponential fit of ~ 0.61 125- The Milky Way profile is an excellent match to the average trend. 



on r 2 5 2 , which varies by two orders of magnitude for the 
galaxies plotted. 

This leads to the surprising result that the mass of 
neutral hydrogen gas depends mostly on the size of the 
stellar disk and that the gas arranges itself somehow into 
a distribution that is self-similar among galaxies. The re- 
sult also implies that except for the region of a galaxy at 
< 0.2 r25, which is highly variable, the total neutral gas 
mass of a disk galaxy at z = can be estimated if 125, 
i.e., the extent of the stellar disk, is known. The variable 
inner regions of galaxies, are, however, quite small, and 
despite either depressions or large excesses of the molec- 
ular gas in these regions, they typically only account for 
a small fraction of the total neutral gas mass (~ 15% on 



average in our sample) . 

We examine this conclusion in Figure SI which shows 
the total gas mass of the galaxies in Figure |3] as mea- 
sured in the respective H I and H2 intensity maps, com- 
pared to the predictions made using Equation [2] (where 
the integration is carried out to 2xr2s). We estimate 
the uncertainty on the total gas masses derived from 
the intensity maps (x-axis) from the calibration uncer- 
tainties of the TH INGS and HERACLES data (5% an d 
20%, respectively. iWalter et aLll2008t ILctov et al.ll2009t) . 
We thus adopt 20% as the overall uncertainty. For the 
masses derived from the fit (y-axis) , we estimate the un- 
certainty from error propagation (from Equation [2]) . We 
take into account the uncertainties on the fit (slope and 
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Fig. 4. — Total (H I+H2) gas mass obtained from the average radial profile fit (Figure[3]and Equation^ versus the mass obtained directly 
from the H I and H2 intensity maps. The flux is computed out to two optical radii V25 for each of the non-interacting galaxies shown in 
Figure [3] The solid line indicates unity and the two dashed lines a factor of two scatter. The plot thus shows that Equation [2] allows us to 
predict the total gas mass of a disk galaxy within a factor of two uncertainty. 



intercept), r 2 5 (estimated from iPaturel et al -111 9 911) and 
Strans (from lLerov et al1l2008D . This results in an uncer- 
tainty estimate of ~ 40% for the mass prediction. Both 
uncertainties are indicated as error bars in Figure 2J 

Figure 2] shows that this equation offers a robust way 
to predict the total gas mass of disk galaxies to within 
a factor of two (indicated by the dashed lines, the solid 
line indicates unity) . The mean ratio of predicted versus 
measured total gas mass is 1.17 (with standard deviation 
0.36), which underlines the good correspondence between 
prediction and data. This close correspondence is rooted 
in the tightly constrained fit in Figure |3l which falls al- 
most directly on the means in the individual radius bins. 

Based on observed metallicity g radients across galaxy 
disks (e.g.. iMoustakas et al.ll2010T ). one might speculate 
about a similarly radially varying CO-to-H2 conversion 
factor. Recent work indeed suggests that systematic vari- 
ations of Xqo with radius (albeit not necessarily driven 
primarily by metallicity) occur in at least some of the 
galaxies studied in this paper: in these cases, Xco is 
observed to increase with increasing radius (Sandstrom 
et al., in prep.). A low value for XCO in galaxy centers 
seems to be a more common feature (Sandstrom et al., 
in prep.). The former effect would lead to shallower H2 
profiles , while the latter might in fact explain the ob- 
served apparent H2 excesses in galaxy centers, possibly 
even giving rise to truly exponential H2 profiles including 
the centers. 

One of the implications of a universal gas profile is that 
there appears to be some mechanism that keeps the re- 
lationship constant in t he face of sta r form a tion in these 
galaxi e s. For example, Bigiel et al. (|2008f ): iLeroy et al.l 
(|2008( k iBigiel et al.l (|2011f > show that the depletion time 



for the molecular gas in disk galaxies is ~ 2 x 10 9 yr. 
The region of Figure 3 where s ? gas > 1 (roughly for 
r < 0.45 x T25) is the region dominated by molecular gas 
and that gas will therefore be exhausted in about 15% of 
a Hubble time. In order to keep the profiles self-similar, 
this implies that either new gas comes from outside the 
disk and falls preferentially in the central regions, an un- 
likely occurrance, or that gas flows through the disk to 
make up the gas lost to star formation, which takes place 
preferentially within r < 0.45 x r25. In either case, it is 
difficult to understand why the total gas profiles would 
be so self-similar. For the infall case, one expects such 
infall to be sporadic, possibly occurring in the form of 
small galaxies merging with the bigger disk galaxies in 
our sample. For the inflow case, Figure [3] encompasses a 
variety of galaxy types and morphologies, and it is hard 
to see how inflow could be so closely regulated to produce 
the tight observed relationship. 

As cosmological simulations of galaxy evolution includ- 
ing gas become ever more refined, it will become neces- 
sary for these simulations to reproduce the universal rela- 
tionship when carried out to the present epoch. Because 
of the limitations in making observations of the atomic 
gas to higher redshift with present day instrumentation, 
it will be difficult in the near future to extend the work 
presented here to normal disk galaxies at significantly 
higher redshifts. It will also be a challenge to extend 
this work to lower metallicity systems, where CO emis- 
sion becomes a n increasingly poor t racer of the molecular 
hydrogen (e.g.. iBolatto et aLll2011[ ). On the other hand, 
r25 can be determined in many galaxies to higher red- 
shifts and Equation 2 could be used to estimate M gas . As 
it becomes possible to measure total neutral gas masses 
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to higher redshifts with ALMA, the JVLA, and Arecibo, 
direct comparisons of the predicted and measured total 
mass of neutral gas can be extended both to higher red- 
shift and to a larger sample of local galaxies providing 
good tests of the universality of the neutral gas profile 
presented in this paper. 

The universal gas profile we have described seems to 
be a fundamental property of normal disk galaxies at 
z = 0. We have, however, only probed galaxies of near 
solar metallicity and none of the galaxies in o ur final sam- 
ple ar e dwarfs. Also. I Young et afl (|2011l ) and lSerra et al.l 
(2012) have shown recently that a surprisingly large frac- 
tion of early type galaxies, i.e., ellipticals and lenticu- 
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lars, contain large amounts of atomic and molecular gas. 
Whether these galaxies obey the same universal gas pro- 
file is still to be determined. 
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